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Abstract 

We apply a gauge invariant formulation of Wilson Renormalization Group (RG) 
to the computation of the Debye and transverse gluon masses in pure gauge SU(N) 
at high temperature. Following the Hard Thermal Loop effective field theory as a 
guideline, we develop an approximation scheme to the exact evolution equations. 

The Debye mass receives sizable corrections compared to the leading order per- 
turbative result, mainly due to the infrared singular behavior in the transverse gluon 
sector. A non-vanishing mass for the transverse gluons is found, which acts as an 
infrared regulator though not efficiently enough as to restore the validity of pertur- 
bation theory. Indeed, discussing the role of higher dimensional operators, we show 
that the gauge coupling for the transverse modes typically flows to non-perturbative 
values unless extremely high temperatures are reached. 

After comparing our results with recent lattice simulations, we comment on the 
possibility of using this formulation of the RG as a tool to construct an effective field 
theory for the non-perturbative, long wavelength, transverse modes. 
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1 . The study of the properties of a non-abelian plasma at high temperature is of great 
interest both in cosmology, concerning the quark-hadron phase transition and the properties 
of the symmetric phase of the Standard Model, and in particle physics, in connection with 
the possibility of creating the quark-gluon plasma in heavy ion collisions. 

From a theoretical point of view, it has been known for a long time that perturbation 
theory is afflicted by severe infrared problems caused by the static transverse gluons, which 
remain massless at any order [Q. As a consequence, the loop expansion ceases to be 
an expansion in powers of the coupling constant and, for instance, infinite loop orders 
contribute to the pressure at the same order 0(g 6 ). 

The Debye (or longitudinal) mass describes the screening of chromo-electrostatic fields 
in the plasma, and its numerical value influences the probability of J/ty and T forma- 
tion, among the main signatures of quark-gluon plasma formation. At leading order in 
perturbation theory it is given by m\ = Ng 2 T 2 /3 for SU(N). At next to leading order 
infrared divergencies show up, and the most one can compute is the coefficient of the term 
g 3 \og(jriL / tut) , ttlt being the transverse (or magnetic) gluon mass, which at this level 
has to be introduced by hand 0. The coefficients of the higher order terms in g are not 
computable in perturbation theory. 

In this letter we will apply the Wilson Renormalization Group (RG) method introduced 
in refs. 0, || to the computation of the Debye mass m L . As we have seen, the value of 
rriL beyond leading order is influenced by the infrared cut-off in the transverse sector, my. 
We will first derive a coupled system of RG flow equations for tul and tjit describing the 
effect of the integration of thermal fluctuations at larger and larger length scales. Then, we 
will include in the system also the running of the coupling constants and of wave function 
renormalizations. 

The RG formulation developed in refs. |3|, |J is based on the introduction of an infrared 
cut-off, A, in the thermal sector of the theory, so that the modes of momentum \k\ 3> A 
are in thermal equilibrium at the temperature T while those with |A;| <C A are frozen at 
T = 0. The RG flow equations then interpolate between the full (renormalized) quantum 
field theory at T = in the A — > oo limit, and the full quantum field theory in thermal 
equilibrium at the temperature T for A — ► 0. Compared to other formulations of the 
Wilson RG in the literature ||, we are then integrating out only the thermal fluctuations, 
all the quantum fluctuations being already included in the initial conditions of the RG 
flow. As discussed in detail in ref. J| the main advantage of introducing the cut-off only 
on the thermal sector is that BRST invariance is preserved. From a computational point 
of view, this allows us to use Slavnov- Taylor (ST) identities as a powerful constraint in 
approximating the exact evolution equations. From a physical point of view, we have a 
tool to derive an effective, gauge invariant, field theory, even for a non-zero value of the 
cut-off A. 

The effective field theory at the scale gT has been developed in refs. and is known as 
Hard Thermal Loop (HTL). As external momenta smaller than gT or next to leading order 
corrections are considered, the magnetic divergencies show up and the HTL resummation 
breaks down. Using the HTL effective theory as a starting point, the RG will guide us 
deeper in the infrared. As we will see, as A approaches tjit the gauge coupling in the 
transverse sector rises to non-perturbative values, then large corrections are to be expected 
to our results for %. On the other hand, most of the renormalization of takes place 
for larger values of A, where the couplings are still small, and the approximation to the RG 
equation is still reliable. 
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2. In this letter we will be mainly interested in the longitudinal (or Debye) and trans- 
verse (or magnetic) masses, defined as the static poles of the full propagator, 

m l,T = n L,r(go = 0, \q\ 2 = -m\ T ) , (1) 

where Hl,t axe obtained from the self-energy IP^ as IT^ = II 00 and Ut = —1/2 IP 1 . As 
shown in ref. |7[], the definitions ([I]) are gauge- independent. The same holds even in presence 
of the infrared cut-off, since it does not break BRST invariance. 

The RG flow of the Debye and magnetic masses is described by the equations H: 



d ___ 2 _ A^U L<T (q) 



A— — mi 



dA L ' T i _i_ 9n L , T (g) 
1 + %T 



(2) 



\q\ 2 =—m? 



L,T 



where the flow equations for Hl,t are obtained by taking the appropriate components of 
that for the self-energy. 

In the case of one particle irreducible vertices the RG flow equations can be obtained 
by the following simple recipe: 1) take the expression for the one loop correction to the 
desired vertex; 2) substitute the tree level propagators and vertices in it with the full, cut- 
off dependent ones; 3) take the derivative with respect to the explicit cut-off dependence 
(i.e. derive only the cut-off function in the propagators). For the self-energy we get Q 

d „ , . if d 4 k 



A—U Ai ^(q) = -- J j—r^Kx, P \(k) [T A ^ pX (q, -q, k, -k) 



(27T) 

+2 T AtPpv (q, -k, k-q) G A>r)S (k - q) T A>1/XS (-q, k,q-k)) + Ghost. (3) 

The kernel, 

K A , pX (k) = p AipX (k)s(k )A5(\k\-A)N b (\k \) , 
contains the Bose- Einstein distribution function, N b (\k \), and the full spectral function 

pA, P x(k) = -ie(k )DiscG AjpX (k) , 

and G AiP \, T A ^ PV , and T A ^ up \, are the full propagator and vertices computed at the scale 
A. 

Up to now no approximation has been performed, and the evolution equations (@) define 
the screening masses non-perturbatively. Of course, in order to solve them, one should know 
the exact propagator and vertices (including their momentum dependence) for any value 
of the cut-off A. This requires managing an infinite system of differential equations for the 
self-energy and all higher order vertices, which is equivalent to solving finite temperature 
QCD. 

As discussed for instance in 0, RG equations can be solved iteratively, reproducing 
the usual loop expansion. Introducing in the RHS the A-independent expressions for the 
propagator and vertices at the loop / and integrating in A yields the results at the l+l 
loop order. In this case, we recover the well known results m^ = mjf = JN/3gT and 



m? = after the first iteration. However it is well known that in high temperature QCD 
the loop expansion ceases to be a sensible approximation already at the second loop, since 
contributions 0(g 3 ^ 2 T) to m^ are generated at each higher order || |2]. Concerning m^, it 
is easy to realize that starting from = 0, no non-vanishing value can be generated at 
any order in the loop expansion. 
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Analogously, the HTL loop expansion at l+l order can be reproduced by introducing 
the A- independent results obtained from the HTL effective theory at the /-loop in the RHS 
of eqs. (fj). After the first iteration, the result of eq. @ for m L discussed by Rebhan in 
ref. is obtained. 

In order to go beyond the HTL effective field theory and closer to the RG framework 
the next logical step is to promote the A- independent HTL propagator and vertices to A- 
dependent ones. Then, the contribution to m^j at the scale A will be given by the RHS of 
eqs. (H) with m L T (and the coupling constant) evaluated at the same scale A. Moreover, 
we will require that the improved vertices and propagator satisfy the same tree level ST 
identities as the HTL ones || . 

The HTL propagator and vertices are complicated functions of the momenta ||, whose 
expressions simplify when the external energy vanishes. Since in this letter we are interested 
in computing the static quantities in ([]]) a great simplification to the flow equations can be 
obtained by rotating to the imaginary time || . The evolution equation (^) then takes the 
form (from now on we omit the A-dependence of the various quantities) 

d 

A— n^(g = 0, \q\) = T^F^z = 2i-KnT;\q\;m L ,m T ) 



dk 

Tr- F iMy( k o'Aq\; m L,m T ) (4) 

Z7T 



with 



1 f d 3 k - 

F^ko] \q\;m L ,m T ) = -- J -^-^AS(\k\ - A)G pX (k) [r^ pA (g, -q, k, -k) 

+2 r MP ^(g, -k, k-q) G vS (k - q) T uX s(-q, k,q-k)} + Ghost . 

Notice that, due to the cancellation between the Matsubara term in the first line and the 
"T = 0" term in the second line, eq. (||) is free from ultraviolet divergencies, as it is 
manifest also in the form (|3|), which contains the Bose-Einstein function. This is of course 
to be expected, since these RG equations describe the effect of thermal fluctuations only. 

The evolution equation for involves F Q0 . We approximate it by using tree level 
propagator and vertices for the non-zero Matsubara modes and in the "T = 0" part, 

F 00 (k ; \q\;m L ,m T ) ~ F 00 (k ; \q\ =0;m L = 0,m T = 0) for n ^ . 

Analogously, in the evolution equation for my, we approximate Fa as 

Fu(k ; \q\;m L ,m T ) ~ F u (k ; \q\ = 0;m L ,m T = 0) for n ^ . 

These approximations will lead to 0(g 2 m 2 L /T 2 ) and 0(g 2 m\./T 2 ) errors, respectively, since 
0(g 2 mL t T/T) terms may only come from the, infrared problematic, zero Matsubara mode. 

The next-to-leading correction in the HTL resummed perturbation theory can be ob- 
tained by using resummed propagators and vertices for the zero mode only 0]. In the same 
spirit, we will employ HTL-inspired propagator and vertices only for n = 0. 

The propagator we need is 

A I Zl 



L 

Z T ( k^kA k^k u 



(5) 

fco=0 
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where, compared to the true HTL propagator [Q] we have A-dependent m LjT and 
Zl^t- The wave function renormalizations are defined, as usual, as Zlt = 1 + 
' 9 ',IT,,, T (g = 0,|gl 2 ) 



W= m L,T 

In the HTL approximation, the three and four gluon vertices, evaluated at zero external 
energy, reduce to their tree level expressions. Here, the presence of a non-zero my forces 
us to modify the trilinear vertex in order to preserve the gauge invariance of the evolution 
equations @. Requiring that the new vertex is related to the propagator (|5|) by the tree- 
level ST identity does not fix it completely The transverse component of the vertex could 
be in principle determined by studying its RG equation. In order to limit the number of 
flow equations, in this letter we will instead fix the transverse part by imposing that the 
improved vertex has the same form as that obtained by Alexanian and Nair M, 



where 
and 



r^ P (g, k, -q - k) \ qo=ko=Q = g[v pvp (q, k, -q-k) + V pup (q, k, -q - k)] (6) 
Vp(?' ~q ~ k) = (q - k) p g^ + (2k + q) il g vp - (2q + k) v g w 



with 



V 0up (q, k, -q - k) 
V ije (q, k, -q - k) 



A 



V^opiq, k, -q-k) = V^q, k,-q — k) — 0, 

-m 2 T A [B ki kj k e + C (g» qj k e + qj q e k { + & q e kj] 

- (q <-> k)} 



k 2 q 2 — (q ky 



B 



q k 



q{q + k ) 

(q + k) 



2 ' 



c 



k(q + k) 
(q + k) 2 



Note that the vertex correction V pup does not appear in the evolution equation for m L but 
only in that for %. Inserting eqs. (|5|), (P), and the tree level expressions for the four-gluon 
vertex and the ghost propagator and vertices in the on-shell evolution equations (0) the 
gauge parameter dependence drops, as can be checked by explicit calculation. 

In the flow equations (0) the vertices are evaluated at external momenta of order A. As 
a first step, the thermal contribution to the renormalization of the coupling constant will 
be neglected, and the dependence on the external momentum will be taken into account 
by using the T = beta function at two loops, 



^ 2 (A) 



11 JV 
24tt 2 



log 



A 



A 



MS, 



17 JV, 
+ 88^ l0g 



2 log 



A 



A' 



MS, 



(7) 



for A > m,L, and g(A) = g(mi) for A < m^. The wave function renormalizations Z^t will 
be fixed at Z^^ = 1- 

The thermal renormalization of ttil/T as a function of the cut-off A can be read from 
Fig.l. The effect of the fluctuations integrated out at scales A/T ^ O(10) is negligible 
due to Boltzmann suppression. Around A/T = 0(1) the HTL contribution is rapidly built 
up. Stopping the running at this scale basically reproduces the HTL effective field theory 
results. Going to lower energy scales the infrared divergencies in the magnetic sector start 
showing up and, depending on the value of tut, their effect may eventually overcome the 
leading perturbative result. This is clearly seen from the upper line, where we have set 
mx = 0. 
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Figure 1: Running of m^/T vs. A/T in SU(3) for T = 10 4 Aj^. The solid line is obtanied 
fixing m T = 0, and the dash-dotted line for m T = g 2 {T)T. The dashed line corresponds to 
the solution of the system of RG equations for m L and tut- 



In next-to-leading order perturbative computations the generation of a 0(g 2 T) magnetic 
mass by non-perturbative effects is usually invoked as a mechanism to regulate the infrared 
divergencies coming from the transverse gluon sector. The computation of at next-to- 
leading order gives [p] 



»n = ^g 2 T 2 (l + ^g\o g ^ + Cg + 0(g 2 )) (8) 
3 V Z7r m T J 

where C is a coefficient receiving contributions at any loop order and therefore not com- 
putable in perturbation theory. If we introduce by hand an ad hoc ttit = 0(g 2 T), then 
perturbation theory can tell us the coefficient of the glog(l/g) term, but in any case it 
becomes meaningless beyond this order. 

Fixing % = g 2 (T)T, where g(T) is obtained by computing eq. ([7|) at the thermal 
scale of perturbation theory, = 47rTexp(— j E ) ~ 7.055 T [|1(J, the lowest line in Fig.l 



is obtained. In Fig. 2 m^/T is plotted as a function of the temperature. The lowest line 
corresponds again to tut = g 2 (T)T. Notice that it agrees with the next-to-leading order 
perturbative result (eq. (§) with C — 0) only for very high temperatures T <; 10 10 A^r. 
The next step is to let rrir run, studying the coupled system of RG equations for m L and 
%. We note immediately that, at this level of the approximation, = is a fixed point 
of the RG flow. In fact, in this case, the RHS of the flow equation reduces to the one-loop 
expression which, upon integration, gives my = 0. 

The possibility of a thermal renormalization of the magnetic mass is then related to 
its initial value at T = 0. The structure of the gluon propagator is still subject of intense 
study, and no definite answers are available. However different groups have found evidence 
for the existence of a pole at non-zero momentum of O(Ajj^) [|11|]. In the following we will 
assume that such a pole indeed exists, and will use my = mi = V-^JTs^ with r\ a coefficient 
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Log,, (T/AJ 



Figure 2: Plot of m L /T vs. log 10 (T/A^-) for SU(3). The lowest line is obtained by fixing 
m T = 9 2 (T)T, the "HTL" line represents eq.(f|) computed for g = g(T) and C — 0, and 
the "HTL + NP" (dotted line) is the same equation but for C = 1.3. The solid line is the 
result of the RG system for and uit- 



of 0(1), as initial value for the running. Increasing T, the sensitivity to rj becomes weaker 
and weaker, since changing this parameter is equivalent to rescaling Ajjg, which enters the 
problem only in the argument of the logarithm in eq. (0). 

The result is plotted with the continuos line in Fig. 2. A comparison with the next-to- 
leading expression of eq. (§) (dotted line), computed for g = g(T), allows us to determine 
the non-perturbative coefficient C as 




1.1 for SU(2) 
1.3 for SU(3) 



As a general fact, we observe that the corrections to are larger than those obtained 
by fixing my = g 2 T. This is due to the fact that the values of my given by the RG are 
generally much smaller than g 2 T, so the decoupling of the RG flows takes place deeper in 
the infrared. It can be seen explicitely from the dashed line in Fig. 1. A tjit = 0(g 2 T) 
naturally emerges from gap equations, |L2j, § where the A-dependence of rriT inside the loop 
integral is neglected. Once this dependence is kept, the simple scaling law tut = 0(g 2 T) 
is lost. We can understand it by looking at the approximate form of the flow equation for 
rriT obtained in the limit A 3> tut, 

K 9 —m 2 T = -Kg 2 T^ , (9) 

with K = 55/87T 2 ~ 0.697 in the case of SU(3). We may also neglect the A-dependence 
of g and the contribution of the running for A below as we can verify by comparison 
with the full numerical results. 
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The gap equation approximation corresponds to integrating eq. @ with a A- 
independent rriT = m^gap such that 

™T,gap m 2 T (k = 00) 2 m T,ga P f°° dA 



+ K & r 

1 J rriT , 



T 2 T 2 ' * T J mT , gap A 2 

„ Kg 2H^L (T > mr(A = oo) = ^A^) , 

leading to the non-vanishing solution mT, ga p — Kg 2 T. 

On the other hand, integrating eq. (^) down to A = m-p without further approximations, 
gives the equation 

mtp = m^(A = oo) exp {Kg 2 TjmT^ , (10) 

which always gives mr <C m T,ga P for T r] A.j^. 

From the full RG equations we obtain values for m-r which are in reasonably good 
agreement with the solutions of eq. flTUD and, in the range of temperatures considered, are 
very well approximated by the law (in SU(3)) 



m T = 0.128 g 2 T(l + 2.26 g log g) 



3. In order to assess the reliability of our approximation to the RG flow equations for 
rriL and uit we must consider the thermal running of higher dimensional operators. In 
particular, we will study the running of the gauge coupling constant in order to determine 
the efficiency of the screening masses as infrared regulators. 

In the flow equations for tul and mx only two types of trilinear vertices appear, T 00i 
and Tijk, where k = 1, 2, 3. The RG equations for these two vertices differ remarkably 
in the infrared, since that for T^k receives contribution from the loop in which the three 
circulating gluons are all transverse, of mass tut, whereas in the loops contributing to the 
running of r oi at least one longitudinal gluon appears. As a consequence, it is convenient to 
modify the parameterization of the trilinear vertex in eq. (^) by introducing two different 
coupling constants, gi for T 00i , and for T^. Of course, due to the presence of the 
thermal bath, this is not in contradiction with Lorentz invariance. 

From the RG flow of T^ up we extract those of gL,T, defined as 



3l = q2 _ k 2 ? oo P [q,k,-q-k] 



q =0 ,\q\ 2 = -m 



2 



9T = }^ 2{q 2 -k 2 ) Vup[q ^- q - k] 



1 



(11) 



go=0 ,\q\ 2 =-m? 



The running of the couplings is shown in Fig. 3. Notice the impressive difference, 
induced by thermal effects, between gi and g?. This is due to the fact that graphs with 
only transverse gluons in the loop contribute to the running of but not the ones of gi- 
When nix < A <C T the transverse sector of the theory is in the three-dimensional regime, 
and the coupling constant g? grows as 1/A in the infrared. The non- vanishing leads 
to decoupling for A ^ m T , but it takes place when the coupling g T has already grown to 
non-perturbative values. 
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Figure 3: RG flows for g L e Qt defined as in eq. ( [TT]) for SU(3) and T = 10 4 Aj^-. 

Comparing Figs. 1 and 3 we see that the evolution of the coupling constants affects 
our results for itil and my in a different way. First, notice that in the flow equation for 
rriL only the vertex r oi, and then the safer coupling g^, appears explicitly. The effect of 
gx is felt by only indirectly, through the running of m^. Moreover, the running of 
m L is almost completely decoupled before starts its crazy rise. As a consequence, the 
non-perturbative regime for g? does not affect very much our results for m^, since there is 
no time enough to communicate it to the longitudinal sector. 

On the other hand, since most of the running of ttit takes place when g? is also running 
the inclusion of higher order operators gives larger corrections to ttit/T. We estimate that 
they are typically of 0(1), corresponding to a few percent corrections to m L /T. 



4. Summarizing, we have employed the Wilson Renormalization Group at finite temper- 
ature introduced in refs. || || to study the gluon self-energy in high temperature QCD. In 
order to keep contact with the perturbative and HTL-resummed results, we have developed 
an approximation scheme based on the use of HTL-inspired Ansatze for the propagator and 
vertices appearing in the RG equations. 

We have concentrated on w,l and my, the Debye and magnetic masses, deriving a 
system of differential equations for them, the coupling constants and the wave function 
renormalizations. We obtain large corrections to with respect to the leading order 
result, = ^jN/3gT, mainly due to infrared effects in the transverse sector. 

Our results for in pure gauge SU(2) are in good agreement with the lattice results 
of ref. jn| on the pole mass of the longitudinal gluon propagator in the Landau gauge. For 
instance, for T/Ajjg = 4000 we find rni/T = 1.12 whereas they obtain m^/T = 1.1 — 1.3, 
the leading order perturbative value being m L /T ~ 0.84. Lattice computations of the 
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screening masses have been reported also by Kajantie et al. in ref. [15|, where the long 



distance behaviour of certain gauge-invariant operators is studied. In this case corrections a 
factor two to three larger than ours and those of ref. |13| are generally found. As discussed 
[13|] a possible explanation of this discrepancy could be that the effective thermal 



m 



mass extracted from gauge-invariant composite operators arises from the superposition of 
several decoupled gluons. 

Concerning %, results much smaller than the value g 2 T obtained from gap equations 
are typically found, and this has been ascribed to the scale dependence of m T . 

The inclusion of the running of the coupling constants typically leads to sizable correc- 
tions in the transverse sector, since the relevant coupling rises to non-perturbative values 
in the infrared. This result clearly indicates that the dynamics of the transverse modes of 
wavelength ^ l/mr is highly non-perturbative. 

In perturbation theory, the non-perturbative modes circulate in the loop at any order, 
leading to the infrared problem that we have already discussed. The same statement that 
the HTL give the effective theory at the scale gT is meaningful only at leading order, since 
at higher orders the longer scales come into play. On the other hand, in a RG framework, 
the scale A can be used to keep the long wavelength modes under control, allowing the 
definition of effective theories at larger and larger scales in a clean way. This procedure 
usually runs into problems for gauge theories, since the introduction of the cut-off A breaks 
gauge invariance. In our framework this is not the case, since the cut-off is imposed on the 
thermal, on-shell, sector of the theory only The construction of effective field theories 
along the lines illustrated in ref. |16| for the scalar theory can then be properly carried out 
in this context. 
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